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Abstract 

This article is intended to provide a pedagogical introduction to the supersymmetry method for per- 
forming ensemble-averaging in Gaussian random-matrix theory. The method is illustrated by a detailed 
calculation of the simplest non-trivial physical quantity, namely, the second-order correlation in the den- 
sity of states for two different energies within the spectrum (commonly known as the density-density 
correlator) for a system described by a random Hamiltonian matrix belonging to the Gaussian unitary 
ensemble. 
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I. Introduction 



In many physical problems, especially in nuclear physics and condensed-matter physics, the quantum 
mechanical Hamiltonian, which models the underlying physical complexity, involves a random matrix 
belonging to one of the three standard Gaussian ensembles: the Gaussian unitary ensemble (GUE), 
the Gaussian orthogonal ensemble (GOE) or the Gaussian symplectic ensemble (GSE). The appropriate 
choice is dictated by the symmetries of the theory 

In condensed-matter physics, such problems arise in the study of electron localization phenomena 
in disordered conductors and semi-conductors. In the mesoscopic regime, applications include universal 
conductance fluctuations [||, Aronov-Altshuler-Spivak oscillations Q and persistent currents ||^. These 
are all weak-localization effects. In the strongly localized domain, problems such as long disordered 
quantum wires and the integer quantum Hall effect |Q are amenable to a random-matrix formulation. 
In the ballistic regime, where disorder is negligible, classical chaos in the dynamics of electron scattering 
gives rise to the complexity that can be described within random-matrix theory. Thus, another recent 
important application has been the study of electron transport across microstructures, constructed in the 
shape of classically chaotic billiards, in the presence of a variable external magnetic field. This allows one 
to investigate the continuous crossover between two random-matrix ensembles (GOE GUE) as probed 
by the mean conductance and its correlations for increasing magnetic field Also, connections between 
the Sutherland-Calogero model and the random-matrix ensembles have been established and discussed 
in the literature ^, ^ . A brief review of applications to nuclear physics can be found in Ref. Finally, 
the latest application of the supersymmetry techniques described here has emerged in elementary particle 
physics, where random-matrix theory has been used to study chiral symmetry breaking in QCD p^ . 

One is often interested in computing ensemble averages of physical quantities involving traces of 
products or products of traces of resolvents of the Hamiltonian. Such quantities include the density of 
states and its correlation functions, products of S-matrix elements, as well as the conductance (by virtue 
of being itself a sum over products of S-matrix elements, as can be seen from the Landauer formula, which 
states that G = {e'^/h)trtH with the matrix t denoting the transmission part of the S-matrix) and its 
higher moments and auto-correlations. In this situation, it is convenient to proceed by first constructing 
a generating function for resolvents and expressing the desired quantities as derivatives of the generating 
function with respect to the source. The supersymmetry method [^s], provides a mathematical tool 
for performing the ensemble average of the entire generating function, thereby avoiding the necessity for 
expanding in powers of the random matrix, and subsequently resumming after ensemble averaging (as is 
done in impurity perturbation theory, for example). Therefore, it is a non-perturbative approach which 
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can provide results beyond the domain of validity of expansion techniques. It is also a useful alternative 
to the method of orthogonal polynomials in many cases Q. Indeed, where coupling to external channels 
is involved, such as in the calculation of S-matrix correlations |Q, orthogonal polynomials have not been 
applied and may not be amenable to such problems. 

Furthermore, non-perturbative results are especially useful whenever a small number of external 
channels, M, is coupled to the random system, because the perturbation expansion here proceeds essen- 
tially in powers of 1/M. This situation can typically occur in the problem of compound- nucleus scattering 
with few open decay channels , or in the conductance problem for electronic microstructures connected 
to few electron- mode leads Q . An extreme case is provided by isolated (M — 0) non-dissipative systems 
in dealing with the so-called zero mode. The problem of persistent currents gives one such example 

The foregoing remarks have hopefully illustrated the fact that random-matrix theory is a subject of 
much current interest, and the supersymmetric approach is playing an increasingly important role in it — 
its use having been very rapidly increasing over the past few years. It is difficult, however, for the novice 
to easily grasp the elements of the supersymmetry method simply by studying the topical literature. 
Applications to real physics problems typically involve calculations rendered long and laborious by the 
presence of parameters inducing explicit symmetry breaking, couplings to external degrees of freedom, 
the necessity for complicated source terms, continuous spatial dimensions, and the large dimensionality 
of the underlying matrix spaces required for evaluation of higher-order correlators. Usually also, other 
extraneous ingredients come into play in a given physical problem, such as statistical scattering theory. 
All these complications tend to obscure the essentials of the supersymmetry method. The purpose of the 
present article is to bridge the gap by providing a primer for more realistic calculations. 

To this end, we consider here a quantum mechanical theory whose Hamiltonian is taken to be purely 
a random N x N Hermitian matrix belonging to the Gaussian unitary ensemble (GUE). By definition, a 
random matrix H is said to belong to the GUE if (i) its diagonal elements i/^^ and the real and imaginary 
parts of its off-diagonal elements -ff^^, for < ly, are statistically independent, and (ii) the probability 
P{H)'DH that the system belongs to the volume element 

TV 

VH^W dH^^ ■ l[ d{ReHf,,)d{lmH^,) , (1.1) 

11=1 n<v 

is invariant under every automorphism 

H ^ H' ^U^^HU , (1.2) 

where U is any unitary matrix, i.e. P{H')'DH' = P{H)'DH. It should be noted that the diagonal elements 
H^^ are necessarily real, as H is Hermitian. 
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Because of (i), P{H) is a product of functions each of which depends on a single variable (viz. the 
statistically independent elements H^^^^, Re i/^i/, Im /i < v). Then, as a consequence of the unitary 
invariancc (ii), as originally proven by Porter and Rosenzweig, it follows that the statistically independent 
elements are Gaussian distributed with a common variance. A derivation of this fact may be found in 
Section 2.6 of Ref. Accordingly, the second moments of the matrix elements are given by 

N' 



and we take them to have zero mean i/^^ = 0. The appearance of the matrix dimension N in Eq. (1.3) 
provides a convenient normalization for the common strength A, while the structure of the Kronecker 
deltas reflects the required statistical independence as well as the Hermiticity of H . For obvious reasons, 
we refer to /i, v as the level indices. The overbar represents ensemble averaging; and since we assume a 
Gaussian distribution of Hamiltonian matrices, the mean and second moments (as specified above) serve 
to define it uniquely. 

The ensemble average of any function f{H) of the random matrix H can be represented as 

J{H) = 1:J VHf{H)c^p[-^TTH^], (1.4) 
with the measure DH as given in Eq. (^]^) and the constant A/" chosen to ensure unit normalization. 



Now let us suppose that the function f{H), whose average we are considering in Eq. (|l.4[ ), depends only 
on the eigenvalues of H, viz., the energy levels Ei, E2, ■ ■ ■ , Em- This will happen, for example, if f{H) 
satisfies the symmetry property 

f{U-'HU) = f{H) (1.5) 

for any unitary matrix U £ U(A^). To see this, let us note that any Hermitian matrix H can always be 
diagonalized by an x unitary matrix Uh from the coset space U(A)/[U(1)]^, i.e. 

H - UhHuUh' , (1.6) 

where iJo is the diagonal matrix of eigenvalues Hu — diag(£^i, E2, ■ ■ ■ , E^). Then we have 

f{H) = HUh'HUh) 

= fidia.g{E,,E2,...,EN)) 

= fiE,,E2,...,EN) . (1.7) 



Moreover, because of the decomposition (|1.6|), we can make the change of integration variables in Eq. (1.4) 



from the matrix elements i?^^ to the A^ eigenvalues i?i,i?2, ■ • ■ ,E]\f and A^^ — A^ real angles that are 
needed to parametrize the space of unitary matrices Uh- 
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This allows us to explicitly perform the N'^ — N angular integrations implicit in Eq. (1.4), and 
thereby express it as 

f(Eu...,EN)^^ l[dE,JiEu...,EN)f{E,,...,EN)e^p{-—Y,E^}. (1-8) 

2—1 i—1 

where J{Ei^ . . . , Ej^) denotes the Jacobian for the transformation of integration variables, and 

Vn= I dflNiU) , (1.9) 



being the integral of the Haar measure dftN{U) over the coset space U(A^)/[U(1)]^, expresses the volume 
of this space. If we proceed to write Eq. ( |1.8| ) in the form 
n ^ 

f{Ei,...,EN)^ l[[dE,f{Ei,...,EN)P{Ei,...,EN) , (1.10) 

then we can immediately read off an expression for the joint probability density function for the eigenvalues 
as being given by 

P{E,,...,EN) = ^JiE,,...,EN)cxp[^^Y.^^} ■ (1-11) 

i—l 

To obtain the two-level distribution function, one integrates over all eigenvalues Ei except for Ei and E2'. 

V f ^ N ^ 

V2{EuE2)^^ n{dE,J{Ei,E2,...,EN)e^v[-i^Y.^^] ■ (l'^^) 



2 — 3 2 — 1 



To calculate the Jacobian J |T^, let us begin by considering the differential of Eq. (1.6), viz. 
R = UHuU-\ which yields 

dH = dUHjyU-^ + UdHjyU-^ + UHudU-'^ . (1.13) 

Consequently, 

U^^dHU = U^^dUHu - HYyU^^dU + diJo , (1.14) 

having used the relation dU^^U = —U^^dU. Now let (dH) denote the exterior product of the differentials 
of the independent elements of _ff , i.e. 

AT 



(dH) = /\ dH^, 

N 

= [/\dH^^} A[/\dH^,AdH;,] 

N 

= {2if^^~^^^^{/\dH^^}/\{/\dReHf,,AdlmH,,,y (1.15) 
given a suitable ordering of the matrix elements. From the linearity of the wedge product, it follows that 

iU~^dHU) = q{U, U*){dH) (1.16) 



for some polynomial q(U,U*) of the independent matrix elements of U and their complex conjugates, 
which can be shown to have unit modulus \q{U, U*)\ — 1. 



On the other hand, by expressing U = {ui, U2, . . . , un ) in terms of column vectors Ui, i — 1,2, N 



(which correspond to the orthonormal eigenvectors of H), and then substituting this form into Eq. (1.14), 



one can easily show that the exterior product of the RHS of Eq. (1.14) simplifies to yield 

N 



iU-^dHU)^[ Jl {Ej - Ekf /\dE,} A{U-^dU) . (1.17) 



l<j<k<N 



After equating Eqs. (1.16) and (1.17), and taking into account Eq. (1.15), we can read off the volume 
element 

N 

VH = dnN{U)J{Ei,...,EN)Y[dE^ , (1.18) 

i=l 

with the Jacobian identified with 

J(£;i,...,£;^) = 2-^(^-i)/2 H {E,-Euf (1.19) 

i<i<fe<JV 

and dQ,M{U) denoting the Haar measure on the coset space U(A^)/[U(1)]^. Consequently, the joint 
probability density function assumes the form 



N 

l<j<k<N i=l 



p{e,,...,em)^cn{\) n (^j-^fc)''^^p{-^E^'} ■ (^■^'^^ 



The constant Cn{\), which is equal to Cn{X) = 2-^(^-1)/^-! Vat, can be determined directly by 
explicitly calculating M and Vn ■ On the other hand, it can also be fixed by using the fact that the joint 
probability density function P{Ei, . . . , En) integrated over all the eigenvalues Ei must yield unity. Thus 
here, for the sake of brevity, we shall take account of this fact and simply derive an expression for C7v(A) 
by appealing to Mehta's integral [|l|, which states that, for /3 = 1, 2, 4, 



i^l"-"^ l<j<k<N 

So, for the GUE value of /3 = 2, we obtain 



(1.21) 



^ = n/ dE.^-^^^^1^:^.-^ n (^.-^^)' 



N 



i2nf/^Y[j\. (1.22) 



The GUE is relevant for systems without invariance under time reversal. Since our aim is to give 
a tutorial introduction to the application of the supersymmetry method for random-matrix problems, 
we shall present a detailed step-by-step calculation of the simplest non-trivial physical quantity in this 



theory, namely, the density-density correlator. Before doing so, however, we shall demonstrate that the 
second order correlator of the density of states at two energies Ei , E2 within the spectrum of the theory 
coincides with the two-level distribution function P2{Ei, E2) that we have just discussed. One physical 
system to which our calculation applies directly is that composed of small disordered metallic grains that 
are subjected to an external magnetic field [|l3j . 

We shall begin our exposition by relating the connected part of the density-density correlator 

to the ensemble average of the product of retarded and advanced Green's 



C{Ei,E2) = [p{Ei)piE2 
functions G^{E), which are given as resolvents. We show that, in the large- limit 

1 



C{E^,E2) 



Re 



G+{Ei)G-{E2) 



(1.23) 



Our discussion will then detail the construction of a generating function Z{e)^ depending on source-term 
parameters e = (ei,e2), with respect to whose components twofold differentiation yields the product of 
an advanced and a retarded Green's function. After having demonstrated how to perform the ensemble 
average exactly with the aid of supersymmetry concepts, and how to extract the large-A^ limit using 
saddle-point techniques, we shall arrive at the result: 

1 92 

e=0 



G+{Ei)G-{E2 



4 deide2 



TT \ ' 



— j / trgfc(3iitrgfcQ22exp 



2d 



■trgiQ 



(1.24) 



assuming Ei + E2 ^ Q for simplicity, where lo = Ei — E2, = uj + ie, and with d ~ nX/N denoting the 
mean level spacing at E = 0. This is an integral over a coset manifold of four-dimensional supermatrices 
Q, which we show how to parametrize and evaluate in order to derive an expression for the connected 
part of the density-density correlator that has been previously established in the literature on random 
matrix theory by a variety of methods |jl|, |l^, |2^, viz., 

1 



G{cu) 



M)2 



7r(5(x) 



sin^ X 



(1.25) 



with X = TTOj/d. Since we assume Ei + E2 = 0, we use the notation C{uj) = C{lu/2, — w/2) for the func- 
tional dependence on the energy difference u. 

II. Density of States 



The density of states, normalized to unity, is given by 

N 



(2.1) 
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With this definition, it is clear that 
Using the identity 

SiE)-- 

we obtain the representation 

Np{E) = 



/ + 00 
dEp{E) = 1 



1 

2Td 



1 

i 



E -ie E + ie 



Tr 



1 1 

-Tr 



E-H-ie E-H + ie 



-[G+iE)-G-iE)] 
--lmG+{E) , 

TT 



where 



C^iE) = Tr- 



1 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



E- H±i£ 

We note that G^{E) corresponds to a retarded Green's function, while G~ {E) is an advanced Green's 
function, i.e., if we introduce the time-dependent Green's functions 

'•+°° dE 



it follows that 



^e-^^*G±(i?), 
-OO 27r 



G^{t)=Q for t^Q 



(2.6) 



(2.7) 



Given Eq. ( |l.4| ) for the definition of the ensemble average and Eq. {2A) for the density of states, we 
can write the density-density correlator as 

p{E')p{E") = JvH Tt5{E' -H)TTSiE" -H)e'^/^^"^' "' 

= / n '^^^ ^(^1' ■ • ■ ' / d^NiU) SiE' - E,)5{E" - i?fe)e-^/2^' S.^i ^' 

1 -i— 1 

= W jI[dE,J{E',E",Es,...,EN)e^p{-^[E'' + E"' + J2E^)}, (2.8) 

i—3 i— 3 

where J{Ei, . . . , i?Ar) is the Jacobian for the change of variables from matrix elements of H to eigenvalues 
Ei and orthonormal eigenvectors Hi, as discussed in the foregoing section; dQ^iU) is the Haar measure 
over the space of unitary matrices formed from the column eigenvectors U — (ui, U2, ■ ■ ■ , un ), and Vn 



is the volume of this space. The final line in Eq. (2^) follows from the fact that the Jacobian is totally 
symmetric under all permutations of its arguments. Eq. (2.8) can be recast as 

N 



piE,)piE2)^CNiX) / l[dE,JiE,,E2,...,EN)e^p[-—J2Ef}, 

'' i=3 i=l 



(2.9) 



8 



from which it is evident, given Eq. ( 1.12| ), that the density-density correlator coincides with the two- level 
distribution function p2{Ei, £2). 



We are interested in looking at the connected part of the ensemble-averaged density-density correla- 



tor, 



C{Ei,E2) = p{Ei)piE2) - p{Ei) ■ p{E2) 



Using the representation (2.4), and the fact that for large- iV 



G^{E^)G^{E2) = G±{Ei) ■ G±{E2) , 



(2.10) 



(2.11) 



we obtain 



CiEi,E2) 



1 



Re 



G+{Ei)G-{E2) - G+{Ei) ■ G-{E2 



(2.12) 



in the large- limit. Hence we shall focus our attention on computing the quantity G^ {Ei)G^ {E2) ■ 
Eq. ( ^.11 ) can be seen to follow from the triviality of the large- iV saddle point (as discussed later on) 
when both imaginary parts lie on the same side of the real axis. 



It is convenient to introduce the average energy E = \ {Ei+E2) and the energy difference 
Lu = El — £2- Then we can write 



El = E + iui , E2^ E -iui , 



where we define Cj — 5 (e — ii^) ■ Thus we have 



E:; )^El2 + iu^L, _^ 



and so 



G+(Si)G-(^2) =Tr 



■Tr 



^ £^12 ~ H + iCoL J \ EI2 - H + iCoL y 22 

The symbol I2 is just the two-dimensional unit matrix. 



(2.13) 



(2.14) 



(2.15) 



III. Superalgebra 

The supersymmetry formalism necessitates the use of anti-comniuting (or Grassmann) numbers. Such 
numbers satisfy an anti-commutative product law, ?7i?72 = ~'n2Vi- Consequently, they are nilpotent, 
rjf — — 0. Complex conjugation can also be defined for Grassmann numbers. Our convention for 
complex conjugation is 

{V1V2)* = rfiV2 . V** = -V ■ (3-1) 
We note that the product of two Grassmann numbers is a commuting number. Thus, a Grassmann 
algebra (or superspace) can be constructed by combining anti-commuting numbers with commuting (real 
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or complex) ones. Linear algebra, analysis and topology can all be introduced on superspaces. For 
an exposition or review, one can consult Refs. [T^-p^ Here, we shall restrict ourselves to defining the 
basic quantities and constructions needed for later use, and presenting our conventions. Because of 
their relevance to the path-integral formulation of elementary-particle theories, commuting degrees of 
freedom in a superspace are often called 'bosonic', and anti-commuting ones 'fermionic'. A differentiable 
manifold with both bosonic and fermionic coordinates is known as a supermanifold. The superspaces 
and supermanifolds that we shall encounter in the present application are ones with equal numbers of 
bosonic and fermionic degrees of freedom. This situation is sometimes described as a perfect Z2-grading, 
and symmetries on such spaces are supersymmetries. 



A=[ R ] ^ (3.2) 



Supermatrices have the form 

f3 b 

where a, b are sub-matrices of commuting elements and a, (3 are sub-matrices of anti-commuting elements. 
Let us label the elements of the sub-matrices by the indices = 1,2,..., and the four blocks by indices 
a, a' = 0, 1 such that 

^pp' ~ "'PP ' pp' ~ PP ' PP' ~ '~^PP ' PP' ~ PP ' K'^-'^l 

In this scheme, where elements of A are represented as A'^^, , the indices a, a' govern the grading of the 
matrix elements, such that a + a' even corresponds to commuting elements while a + a' odd to anti- 
commuting ones. We shall refer to a, a' as the 'graded indices'. The interchange of supermatrix elements 
can now be compactly expressed as 

Ap^B^/J = (-l)("+"')(^+''')i?f,^'A^;; . (3.4) 

One commonly calls A^^ and A^^ the 'boson-boson' and 'fermion-fermion' blocks, respectively, while A"^ 
and A^^ are known as the 'boson-fermion' and 'fcrmion-boson' blocks, respectively. 

The graded trace (or supertrace) of a supermatrix A is defined to be 

trg A = tr a — tr 6 , (3-5) 

or, equivalently, 

trgA = ^(-l)"A^;. (3.6) 

p.a 

The graded determinant (or superdeterminant) can be defined from the graded trace according to 

detgA = e"si'^^ . (3.7) 

From this, one can show that 
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Constructed in this way, the graded determinant has all the usual properties, such as 
detgA~^ = (detgj4)~^, detg(^-B) = detg A- detg-B, etc. Also, for a supermatrix diagonal in the graded 



indices, A^p' = we have 



. det^(o) 



Consequently, detgl = 1, but does not exist for the zero supermatrix. 

Supcrmatriccs act on supcrvoctors ip, whose components wo denote by (p^, so that a = 0,1 corre- 
sponds to commuting and anti-commuting elements, respectively. Thus, if (p and sxe supervectors, 
then 

= (-i)««v«,v« . (3.10) 

The transpose (sometimes referred to as supertranspose) of a supermatrix A is defined by the re- 
quirement that ip'^A'^(p = {Ail))'^(p for any supervectors (p, tp. It follows that 

(3.11) 

This property can be equivalently expressed as 

{A^Y^, = (-1)"(^+«'U«;« . (3.12) 

We see that double transpose is not equivalent to the identity map. Instead, we have 

A** = kAk, ^** = k^, 

ATT ^ j^^^^ ^tt = A, ^^-^^f 

where k^p' = {—l)'^5°'"'6pp> and the Hermitian conjugate is defined as the composition of complex con- 
jugation and supertranspose, viz., At = A*T. 







' aT /3T \ 









IV. Supersymmetry Formalism 
A. Basic Ideas 

We look for a generating function for traces of resolvents {E^ — H)~^, where powers and products are 
obtained by differentiation with respect to a source. As a simple example, let us consider 

, det(£; - + el) ,^ 



Then 

dZ{e) 



= T^-Fi Z{e) . (4.2) 



de E-H + e 
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All powers and products can be generated if one promotes E — H + el, E — H to block-diagonal matrices 
by taking s — diag(ei, £2, • • ■ , £m), E ~ diag(i?i, i?2, • • ■ , Em), and differentiating with respect to the 



somce components £k, k — 1,2, ... , M. The ratio of determinants appearing in Eq. (4.1) ensures that the 
generating function is properly normalized, Z{0) = 1. It follows immediately that if we can ensemble- 
average the generating function, then differentiation with respect to the source will automatically generate 
ensemble-averaged resolvent products, such as those required by Eq. ( 2.15| ). 



For a Gaussian-distributed random Hamiltonian H , this ensemble average is easy to perform if its 
appearance in Z{s) is purely exponential. Indeed, the central element of the present approach is the fact 
that if _ff is an A'^ X random matrix belonging to any Gaussian distribution with zero mean, then 



e-^TrHA^g-i(TrHA)^ (4.3) 

for any fixed N x N matrix A. This result can be easily demonstrated as follows: Since H — 0, expanding 
the exponent yields 



n=0 



{2n)\ 

i~iy 



= Ekir-^pai'-Hi(Tri?A)2j , (4.4) 

71=0 ^ 

where we have appealed to Wick's theorem, which expresses the fact that for any Gaussian distribution, 
all moments factorize into second moments, i.e., the fully connected parts of all higher moments vanish. 
Thus we see that Afpnirin) denotes the number of ways to divide 2n distinct objects into n pairs, 

= (2n-l)!! ^ Ml. (4.5) 



Eq. (4.3) follows immediately. Such exponential dependence on H can be achieved by expressing the 



determinants in the denominator and numerator of Eq. (4.1) for the generating function as Gaussian 
integrals over commuting and anti-commuting variables, respectively. We shall now proceed to discuss 
this construction. 

Let Aq be a non-singular N x N matrix whose Hermitian part is positive definite. Then one can 
express the inverse determinant of Aq as a Gaussian integral in the holomorphic representation [ p!9| , given 

by 

1 f ^ dz^ ^_,t^„, ^ 



det^o J 27ri 

/c— 1 



where z denotes a complex N-dimensional vector z — (zi, Z2, . . . , zn)"^ and is its Hermitian conjugate. 
There exists an analogous Gaussian-integral form for the determinant of any N x N matrix Ai in terms 
of Grassmann variables . It is given by 

detAi = / '[[{27T)dr]*edT]i 6-"^^^" , (4.7) 
12 



(4.8) 



assuming that the Grassmann integrals are normahzcd according to 

(We shall find this choice useful later on.) As a consistency condition, Grassmann integrals necessarily 
satisfy 

J dr]e = J dr]*i = . (4.9) 

The simplest way of seeing why these identities must hold is to note that the value of an integral (being 
a measure) must be a c-number (i.e. an even element of the Grassmann algebra); however, the expres- 



sions in Eq. (4.9) are of first order in the anti-commuting differentials drji^di]^ (i.e. odd elements of the 
Grassmann algebra). Thus, zero remains as the only consistent value. In the exponent above, ry denotes 
an A'^-dimensional vector of Grassmann components 77 = (771,771, . . . ,rii\f)"^. We remark that the multi- 



dimensional Gaussian integrals (4.6) and (4.7) form the cornerstones of the path-integral formulations 
for systems of bosons and fermions, respectively. Derivations of these identities appear in most modern 
textbooks on quantum field theory and many-body physics, such as Ref. uyi. Nevertheless, for the sake 



of completeness, we present derivations of Eqs. (4.6) and (4.7) in Appendix A 



We can combine the matrices Ao,Ai into a block-diagonal matrix A ~ diag(Ao, v4i). This matrix 
can formally be considered to be a supermatrix (with vanishing anti-commuting blocks) . Then 

dot Af) 



detg A = 



det^i ' 



and we can write 



» AT N 

detg"^ ^ = / n ^ d'^kdzl Y\ drjidi] 

•'7 1 /> 1 



(4.10) 



(4.11) 



fe=i e=i 



where tp is the 2iV-dimensional supervector ip = (2,77) . It is straightforward to show that the represen- 



tation (4.11) generalizes to supermatrices with non- vanishing anti-commuting blocks. 



B. Generating Function 

Let us consider a generating function given by 

Z{e) = T)ctg-^[D + J{e)] 

= exp{~Trgln[i5 + J(e)]} , (4.12) 

where the symbols 'Detg' and 'Trg' denote the graded determinant and graded trace as defined in the 
supersymmetric formalism of Ref. |l^ and the foregoing discussion. Here, the inverse propagator D has 
been extended to a 4N x 4N supermatrix, 

D ^ {E - H)U + iCjL , (4.13) 
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where L^p = {—l)P'^^Spp>S°'°' is the diagonal supermatrix that distinguishes between advanced and re- 

D+ = E -H + iOj , D^=E-H-iOj (4.14) 



tarded parts of D. We have = (5""' [diag(i:i+, D\)]pp, , where 



so that D- = d\, and = 1,2, with p ^ 1 referring to the retarded block and p = 2 to the advanced 
block. The indices a, a' = 0, 1 determine the grading (with a ^ for the commuting (bosonic) compo- 
nents and a ~ 1 for the anti-commuting (fermionic) components). Our notational conventions in writing 
traces and determinants are as follows: Those beginning with lower-case letters (e.g. trg, detg) imply 
summation only over graded indices a, a' and/or block indices p,p' , while those beginning with an upper- 
case letter (e.g. Trg, Tr, Detg) indicate that there is also (or perhaps solely) a summation over the level 
indices /i, i^. 

The source supermatrix J(e) depends on a set of parameters Em labelled by some (multi)-index to, 
and is taken to have the general form 

J(e) = e„,M„, (4.15) 

for some set of 4iV x AN supermatrices Mm(/i, i^). It then follows that 
92 



-Zie) 



- Trg Af„Z?-i- TrgMnD-' + Tvg M^D-^ MnD~' . (4.16) 

e=0 



For the problem at hand, we shall make the choice of source matrix (here independent of the level indices 
pi, I/) 

Jp7'(ei,£2)= Yl Smlpp'Mk^'^', (4.17) 

m=l,2 

where 

Ipp'{m) = Spp^Spm (4-18) 

is a projector onto the p-block with p = m and fc"" = (— 1)"(5"" . In this case, we have J = emM„i with 
the identification 

Mm^ Iim)(g)k , me {1,2}. (4.19) 



With Mm given by Eq. (4.1£), the second term on the RHS of Eq. ( 1.16| ) vanishes if we differentiate with 



respect to £m = £1 and e„ — £2, and we see that 



TrtrM - / 2Tr(i:'+) ^ for to^I , . 

LrgMmV -|2Tr(i?_)-i for to = 2, ^^■^^> 



noting that trg k — 2 here. Therefore, 



1 d 



2 



We shall sometimes employ the obvious notation G{E^) = G^{E). 
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(4.21) 

e=0 



C. Ensemble Average 



As discussed in the previous section, the generating function can be expressed as a Gaussian superintegral, 

Z{e) = J VipVTp e'^'^'^'-^'^ , (4.22) 



where 



£i(^;J)=^(^^,[(i? + JVp, (4.23) 



and the measure T>(fT>(f denotes 

V^V^= n ^^d^;i^i)d^^i^i) . (4.24) 

We employ the notation 

(F,G)^^F(m)G(/.), (4.25) 
and </3p (a*) is a four-component supervector field. The adjoint supervector is defined by Ip — ip''s with 

s^p' = s;6""'Spp' , s; = (-l)(i-")(i+p) . (4.26) 
The presence of the supermatrix s in the definition of the adjoint ensures the convergence not only of 



Eq. ( 4.22 ) but also of the final integral representation of the ensemble-averaged generating function as 
a supermatrix non-linear a-model and ensures the correct combination of compact and non-compact 
symmetries therein p^ . By construction, we have Z{0) = 1. 

The ensemble-averaged generating function reads 

Z{^^ fv^V^exp[-iJ2{KAHp)^)} exp{z ^(^^, [(£;! + z^L + J)^]P} , (4.27) 

p. a p,a. 

and we can write 

exp{-tY,{K:iH^)^)} =cxp{~iTrHs} , (4.28) 
where we have introduced the ordinary N x N matrix in the level indices ^, i/, 

We can use Eq. (^^) to obtain the result 

exp|-iTri/5| = exp|-i(Tri/5)2^| 

= exp{-^trg52} , (4.30) 
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where 

^p7'^E^p(^)^p'(^) (4.31) 
is a supermatrix, and it is straightforward to show that 

Tr 52 = y v)S{v, n) = trg . (4.32) 

The ensemble-averaged generating function can now be expressed as 

Z(i) = j V^pVTpe'^"^'^'-^^ , (4.33) 



where 

2iV 



z£2(5'; J) = trg + itrg(£;i + «ciL + J)S'. (4.34) 



D. Hubbard-Stratonovich Transformation 

Ensemble averaging of the generating function Z{e) has introduced a term quartic in the supervector ip 
into the exponent. Consequently, the supervector integration can no longer be performed exactly. How- 
ever, considerable simplification can be achieved through a Hubbard-Stratonovich transformation which 
serves to eliminate the quartic interaction in favour of new (composite) degrees of freedom comprised of 
4x4 supermatrices cr, which couple to the dyadic form ip{fi)Tp{fj,). It turns out that the expectation of 
(an appropriate graded trace over) a is proportional to the average density of states, and hence non-zero. 
We make this explicit in Appendix C. Thus, it is analogous to an order-parameter field signalling the 
occurrence of spontaneous symmetry breaking, the density of states being the actual order parameter. 
Moreover, the cr supermatrices can easily be decomposed into massless (Goldstone) modes and massive 
modes which conveniently decouple in the large- TV limit, and leave us with a theory of interacting Gold- 
stone modes. These are the degrees of freedom in terms of which all quantities in the final effective 
theory can be expressed, and the Hubbard-Stratonovich transformation serves to extract them from the 
original ones (p{fJ.). Hence, another advantage of this procedure is that the number of degrees of freedom 
is reduced by a factor N, since the cr modes are independent of the level indices /z, v. 

To implement the Hubbard-Stratonovich transformation, we follow Ref. ^ in defining 

iC* (S) = trg -Co ArgLS (4.35) 



and 



W{a,S)^~C* [S-'-^a] . (4.36) 
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where uj-c = Reel;, and correspondingly cjj = ImtD. Then we can write 

'Z{e) = J Va e'^'^"--^'^ , (4.37) 

where 

g^£3(-;J) = f p^p^g^lV(<T,S)+^£.(S;J) ^ (4 33) 



assuming that 

J Va e'^(^^^) =Jva e^^^^^'^) = 1 , (4.39) 



i.e. shifts are allowed. The choice of W{a, S) given in Eq. ( f.36| ) serves to eliminate the quartic dependence 



on if in the exponent of the generating function implied in C2{S; J), since 

iW{a,S)+i£*{S) = -^trga^ _ !^ trg crL - iA trg ct5 . (4.40) 
2 A 

The remaining quadratic dependence on (p is amenable to exact integration, which yields 
g,£3KJ) = exp{-^trga2-^trgai} y P^P^ exp{z ^(^^, [(i?l - (iiL + J - Aa)^]P } . 



= exp|-ytrgcr^- ^^-^trgcri-7Vtrgln[£;i-a;iL + J-A<t]| . (4.41) 
At this stage, it is advantageous to perform the shift of integration matrix cr i— > ct' defined by 

- Act' = -^iL + J ~ Xa (4.42) 
in order to move the source matrix and a)i out of the logarithm. Noting that 

trgcr^ = trgcr'^ + -^ [-2Awitrg(T'L + 2Atrgcr'J-2(DitrgLJ + trgj2] , 

1 

trgcri = trgcr'L+ -trgLJ , (4.43) 
A 

we are left with 

Z{e) = / Va exp\ trga'^ - ^-^ trgaL- ^ trg {a + iX~^ojL) J- trg - iV trgln[£;i - Act] | , 
J ^2 A A zA J 

(4.44) 

where we have dropped the primes on the a', and we observe that lu — lu^ + iuji has been reconstituted. 

E. Integration Supermanifold 

In the absence of symmetry breaking {uj = J = 0), Ci{ip; J) is invariant under linear transformations 
if I— > Tip which preserve the bilinear form Tpip. These 4x4 supermatrices must satisfy 

= sTtfi , (4.45) 
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and constitute a supergroup with compact — SU(2) — and non-compact — SU(1,1) — bosonic subgroups. 
The transformation property induced on S is given by 

S TST'^ . (4.46) 

It is also useful to note that S't = sSs. The generating function expressed in terms of the supermatrices 



(T, as given by Eq. (4.37), should also possess invariance under these transformations in the absence of 
symmetry breaking. The structure of W{a, S) then implies that a should also transform like 5, i.e. 
(7 I— > TaT^^. Hence, the domain of integration over the supermatrices a should span a space invariant 
under a TaT^^ ||2^, Such a space is clearly furnished by matrices of the form <t = T^^P^T, where 
the Pd are diagonal matrices with real boson-boson and imaginary fermion-fermion elements. However, 
this choice still does not lead to a convergent integral ( 4. 391 ) over a. To remedy this, we can add an 
imaginary part to Pd according to 

a = T-'^{Pd + rl-iAL)T , (4.47) 

where A is an arbitrary positive constant. It is advantageous to choose A such that a saddle point of 
£3(0"; 0) lies within the cr-integration manifold. This criterion will fix A uniquely. We have also included 
a real term proportional to the unit matrix because, as we shall see, though not required for convergence, 
its presence is also necessary to allow the saddle point to lie within the a-integration manifold. 



It is convenient to re-express Eq. (4.47) in a slightly different form. We perform a coset decomposition 
of the supergroup elements T = RTq where the R span the subgroup that commutes with L. Then we 
have 

a = To\SP + rl - iAL)Ta , (4.48) 

where SP = Rr^P^R. The supermatrices 6P are block-diagonal in the indices {p,p'), and have real 
boson-boson and imaginary fermion-fermion eigenvalues. Now, the generators of the coset elements Tq 
have vanishing entries in both of the diagonal p-blocks. Thus, we can write 

To - exp ( 7 ) , (4.49) 



where, from Eq. (4.45), we deduce the relationship w = —kw'^ between these two 2x2 supermatrices. It 



follows immediately that LTqL = Tq^. 

If we write a — T^^PTq, then the supermatrices P represent massive modes which decouple in the 
large- iV limit. To achieve this decoupling, we shall observe that one can choose P = (Tq + 5P, where 
(To is the unique diagonal saddle point that lies within the integration manifold, in which case 6P are 
interpreted as the massive fluctuations around it. Thus we have a = ctg + T^^SPT with ctg = T^^aoT . 
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The massive fluctuations can be integrated out in the hniit N ^ oo, which leads to an expression for 



Z{e) identical with Eq. (4.44) except for the replacement (t — s- (Tg everywhere. 



F. Convergence 



Let us now consider the convergence properties of the cr-integral in Eq.(4.39), setting S = 0. We have 



^Wia,0) = -ytrg(.+ fL)^ 



N 

T 



trg(Pd + rl - lALf + ^ trg PdTLT-^ + ^ tvg LTLT-A . (4.50) 

A A J 



Integration over the elements of Pd is clearly convergent: The boson-boson elements are real, while the 
fact that the fermion-fermion elements are imaginary compensates for the minus sign coming from the 
graded trace. Thus trgP| is positive definite, and the shift by the constant term rl — iAL does not affect 
convergence properties. Integration over the anti-commuting blocks of T cannot lead to any divergence, 
and neither can integration over the fermion-fermion block of T since it spans the group SU(2), the 
parameters of which are all bounded (being real angles) . On the other hand, the boson-boson block of T 
spans the non-compact group SU(1,1), which has one unbounded real parameter. However, the final term 
in Eq. ( 4.5C| ) furnishes the damping which ensures convergence, provided Alu^-/X is positive. The middle 



term in Eq. ( 4.50| ) does not cause any problems. In the boson-boson block, it is essentially imaginary 



and so gives an oscillating contribution to the integral. On the other hand, the imaginary elements of 
the fermion-fermion block combine with the imaginary prefactor to give rise to a potentially divergent 
contribution. However, since the fermion-fermion block of T^^LT is bounded, this term cannot compete 
with the damping due to the fermion-fermion part of trgP^. To this end, we note that 



A 



trgp2 + ^ trg PdTLT-' = trg Pd + '-^TLP-' ~ . (4.51) 



A 



2 



This would not be so if the fermion-fermion block of T were non-compact. Furthermore, a non-compact 



fermion-fermion block, which would result from taking s = L in Eq. (4.26), would also cause the third 



term of Eq. (4.5C ) to give rise to a divergent integral, owing to the relative minus sign between the boson- 
boson and fermion-fermion blocks in the graded trace. On the other hand, it is clear that the boson-boson 
block of T cannot be arranged to be compact as well as this would necessitate choosing s — 1, which in 



turn would lead to a manifestly divergent integral over tp in Eq. (4.2S). Including non-zero S does not 
alter the arguments presented above. The general conclusions which can be drawn are that the non-zero 
boson-boson elements of Pd must be real and the fermion-fermion elements imaginary, while the boson- 
boson block of T must be non-compact and the fermion-fermion block compact. Convergence cannot be 
achieved in any other way. 
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G. Saddle-Point Equation 



To determine the saddle points of the exponent of Z{e), as it appears in Eq. ( 4.44 ), we neglect uj (as it 
is assumed to be of order 0{N^^)) and the source matrix J (which can be considered as infinitesimal). 
The saddle-point equation then reads 



da 

or equivalently, 



— [i trg + trg ln(£;i - Act)] = , (4.52) 



Xa^ -E(T + X = . (4.53) 



The unique diagonal solution of this equation assuming a form consistent with Eq. (4.47), namely 
ffQ = rl — iAL, corresponds to 



E . / / E^^ 



■ 2A' ^^V'^UJ ' ^'-''^ 
and we note the consistency condition \E\ < 2A. As a function of E, A{E) reproduces Wigner's semicircle 

law; it is in fact proprtional to the average density of states. We can use the normalization condition 



(2.2) to deduce that p{E) — A{E)/{tiX). It follows immediately that the mean level spacing at £' = is 
given by d = ttX/N. 



H. Decoupling of Massive Modes 

If we scale the SP such that a = ctq + N^'^^'^T^^SP'T, then we have 

trg = tvgcrl + N-HvgiSP')'^ + 2N^'^/'^ trgaoSP' , (4.55) 

and 

trgln(S - Act) = trgln(S - Actg) + trgln(l - N-^^^aoSP') (4.56) 
by the saddle-point equation ctq = X{E — Actq)^^. Thus, 

- y trg 0-2 _ AT trg ln(£; - Act) (4.57) 
= -f trgCT^ - iVtrgln(£; - Actg) - ^ trg(<5P')' + ^ trgCTo<5P'CTo<5P' + 0{N-'/^) 
= trgCT^-7Vtrgln(£;- Actg) trg [l-t-i(A2-r2)-firAL] (SP')^ + OiN-^^^) . 



Also, 



trgCTi = trgaGL + N-'^/^trg6P'TLT~^ 

trgCTGi , (4.58) 



N-*c 
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since, given that NCj ^ 0{N^), the latter term wiU produce only non-leading contributions. Similarly, 

trgffj = ivgcjGJ + N-^/'^iYgSP'TJT-^ 

~ trgacJ, (4.59) 

since source differentiations on the latter term will produce contributions suppressed by 0{N^^/'^) com- 
pared with those coming from the first term. Hence, we obtain 

iC^ia;.]) - -itrg[l + i(A2-r2) + irALl (5P')^ + *'C3(ctg; J) . (4.60) 

— >oo 

We see that the variables gq and 5P' decouple in the large- iV limit. 

We note that the space spanned by the supermatrices cr is a linear space. Thus, the integration 
measure can be taken simply to be 

= n '^^p' ' (4.61) 

in which case, complex conjugates of matrix elements, such as (<Jppi )*, should be expressed in terms of 
the original variables, viz. Cp/^- Let us now recall that we have the decomposition a = T^^^PTq, where 
P — 5P + rl — iAL. If we take the independent elements of P and the coset supermatrices Tq as our 
integration variables, then the measure on the a becomes 

Va = I{P)VPVfi{Ta) , (4.62) 

where V^^Tq) is the invariant measure on the To-coset manifold, 

VP = VPi-VP2 , VPp = -i n ^^p"' ' (4-^^) 

for p — 1,2, and 



HP) 



(A?-A^)(A}-Ai) 



(4.64) 



.(A?-Ai)(Ai-AO)_ 

where A^, a = 0, 1, are the two eigenvalues of Pp, for p = 1,2, respectively. Now cr, and hence P, are 
diagonalized by 

cr = T-\Pd +rl- iAL)T . (4.65) 

Thus we see that 

diag(A?, AU", A^) = iV-i/2p, (4 gg) 

after scaling, i.e. 

^r + i{-l)PA + 0{N-^/^) . (4.67) 

Consequently, I{P) = 1 + 0{N^^^^). The massive fluctuations SP{,SP2 can now be easily integrated out 
in the large- iV limit, to yield 

Z(i) ~ K1K2 [vfi{Ta)e'^''^^^'^\ (4.68) 

N^oo J 
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Kp = / I?Ppe-5=p*'-s(^^p)' , (4.69) 



where 



with Cp = 1 + i(A2 - 7-2) - i(-l)PrA for p = 1, 2. Clearly, ii'i = = 1- We show this explicitly in 
Appendix B. Therefore, we see that the large- limit is obtained by setting the massive fluctuations 5P 
to zero and dropping the P-integration, which effects the substitution a —t ctg, as mentioned previously. 
In other words, we treat the massive P-sector at tree level. 



We note that trgo-Q — 0, and by the saddle-point equation ( 4.53| ) 



trgln(£;i - Actg) = - trgln(A"VG) = . (4.70) 

Having set = 0, we see that the correct choice of diagonal saddle-point matrix is ctq = ~iL. So we 
write ctq = —iQ where Q = T^^LT. The supermatrices Q identically span the coset space of the Tq; one 
can replace the coset elements Tq by the group elements T in the expression for Q, as it is clear that only 
group elements from different cosets give rise to distinct Q. We now arrive at the large- A form of the 
ensemble-averaged generating function as a zero-dimensional supermatrix non-linear a-model. 



Z{e) ^ J T^Q e'^-^C^^e'^— =^5'^) , (4.71) 

where 

i£cs{Q) = -^trgiQ, 
iN 

iCsourcciQ; J) - —trg{Q~X-^uL)J, (4.72) 

and it is convenient to write DQ = 'D^(Tq). In the present application, we have trg = 0. Also, we can 
neglect <I> in the source Lagrangian because it will give rise to contributions that are non-leading in A"^. 
Finally, we let oj^ = oj + ie to obtain —Nuj/X = iNuj^ /{2X) = ittuj^ /{2d). Hence we take 

iCesiQ) = ~~2d~^^^^^ ^^(ti-gQii - trg(522) , 

iN iN \ - 

iCsourcciQ;J) = — trgQJ = — ^ EptrgkQpp . (4.73) 

p=l,2 



Therefore, 



1 



e=0 



1 



VQ trg fcQii trg kQ22 exp <^ — — (trg Qn - trg Q22) ^ • (4.74) 



(2A)2 J ^.^.^^..v.^,.^.....^ ^ 2rf 
It is also useful to note that from Eq. ( [4.49 ) it follows immediately that trgQ22 — — trgQu. 
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V. Superintegration 



The measure VQ is the Haar measure on the coset space of the supermatrices Q. If we set e = 0, then 
since by construction Z{0) — 1 and /Csourco(Q; 0) = 0, wc have 

/I 
2?Qexp{^trgLQ} (5.1) 

for any w"*" ^ 0. By now taking the hmit — > 0, we see that VQ should be normahzed to unit coset 
volume. 



A. Coset Parametrization 



Since 





w \ 




w ' 




\^ w 


j 


\^ w 


J 


-i 



WW 
WW I ' 



Eq. (4.49) can be expanded to yield 



V WW 



Let us set 



ti2 = iw ■ 



sm \/ww 



, i2i — iw ■ 



sm V WW 



'WW 



'WW 



so that ti2, t2i are 2x2 graded matrices satisfying 
Then we obtain the representation 

^0 = 



VI + *12t21 i 
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^21 VI + t21<12 

for the general coset element Tq, which leads to 



= 



1 + 2ti2i21 2ti2\/l + ^21^12 

-2i2lVl+il2t21 -(l + 2i21<12) 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



Now let us write 

ti2 = U~^flV , i21 = v~^Jlu , 
where fx,]! are diagonal matrices. This serves to diagonalize ti2t2i and i2i^i2 5 



(5.8) 



ti2t21=U ^{n^J.)u, t2lti2=V 



(5.9) 
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Such a diagonalization must exist because (i) ti2t2i is Hermitian and hence diagonahzable, while (ii) 
^21^12 is not Hermitian but has the same eigenvalues as ii2i2i and so can also be diagonalized. Now, from 



inserting the representations (5.9) into the relations 

(tl2i2l)^ = ^12*21 , (<2lil2)^ = kt2lti2k , (5.10) 

respectively, and noting that the elements of /i/I = flu are real, it follows that one can choose 

u-^^v), y-^^kv^k. (5.11) 



Then, for the representation (5^) to satisfy Eq. ( |5.5| ), we must have fi' = kfi, so that fi, ji can be regarded 



as diagonal special cases of the general matrices ii2,t2i, respectively. For Tq, Eq. (5_^) leads us to the 
form 

To = U-^nU , U=( I M , (5.12) 



where is block-diagonal in the graded indices a, a'. 



> ^ VI + Ai^i 



Thus, let us consider the (bosonic) subgroup of the supergroup of the T that is block-diagonal in 
the graded indices a, a'. Its elements assume the form (rd)^' = <5""'T^p?, and we find it convenient to 



introduce the notation Td = diag(rb,Tf). The relation (4.45) implies that 



T^^ = LTIl, T^^^T}, (5.14) 

where L = diag(+l, —1) in the p-block indices. We see that the subgroup which operates in the boson- 
boson block is Tb e U(l,l), while in the fermion-fermion block, Tf G U(2). However, to generate distinct 
coset elements A = T^^^LTd, we should restrict the subgroup elements to range only over the cosets 
Tb e SU(1,1)/U(1) (a hyperboloid), and Tf e SU(2)/U(1) (the unit sphere). 

These coset spaces can each be parametrized by two real 'angles', according to 

( cosh i^o -Ce*^" sinh i^o \ rj. _ ( cos i^i <e*"*i sin i^i \ 

Ce-'"^" sinh 16*0 coshi6lo ^ \ iQ* e-"*'^ sui\ei cos 16*1 ^^^^> 

where Q is an arbitrary phase, |^| = 1. Let us now combine these parameters into 2x2 (super)-matrices 
given by 

(5.16) 



Then we can write 



00 \ f)_ 
01 j ' ^ ~ I 01 



, 0\f cos^e iCsini^ \ / e"*^ \ 
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where the displayed block structure now pertains to the p-block indices. 



From Eq. (5.13), we can make the identifications 

fi = i(e'^sin^e , iCe-'^ sin ^0 (5.18) 



for the diagonal matrices fi, fi appearing in Eq. (5.8). On introducing the 'eigenvalues 



Aq = cosh 00 , 1 < Aq < oo , 

Ai = cose*! , -l<Ai<+l, (5.19) 



we see that 



Consequently, we have 



Similarly, we obtain 



A^oMo — A'oMo — ^('^0 — 1) , 

HiJI^ = JI^ni = i(Ai - 1) . (5.20) 

giQtrg LQ _ g2ia trg (12*21 _ g2M(Ao-Ai) 21) 

trgfcQii = trgAiifcu^^ , 

trgfcQ22 = -~tTg\vkv-\ (5.22) 

where A = diag(Ao, Ai). We also note that the ensuing representation for Q is then Q — U^^AU with 

X_ { e'^ \ { cose iCsin^ \ / e""^ \ ^on 
1 )[ -zCsine -COS0 )[ 1 J ■ 

The coset supermatrices Tq have two complex Grassmannian degrees of freedom (say rj and p) , one 
of which must be contained in m, and the other in v. These matrices have no bosonic degrees of freedom 
because these are all already accounted for in the bosonic subspace spanned by the Tj. The relations 

= and = kv'^k then imply that we may take 

u - cxdI" - { ^-h*-^ 



~ip* \ _ f 1 + \p*p 



2' 

It follows that we also have the Hermiticity properties = kuk and = v. Thus, 

trgAufcu^i = (Ao + Ai) + 2(Ai - Ao)77*?7 , 

trgAwfcw"^ = (Ao + Ai) - 2(Ai - Ao)p*p . (5.25) 
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B. The Measure 



To calculate the measure 'DQ in terms of specific manifold coordinates, we first note that the invariant 
line element on the Q-manifold is given by 



1 



^tTgdQdQ = ^tig[L,TodTQ^]^ 



1 

16 

i trg [d{h2t2i)d{ti2t2i) - d(ti2VTThIt^)d{t2iVTThIt^) + (1 ^ 2)] . 



(5.26) 



Let us now introduce ri2, T21 by 



tl2 = 2(1 - T12T21) V12 , t21 = 2(1 - T21T12) V21 , 



to obtain the 'rational' parametrization of Tq: 

_f (1 + n2T2l)(l - Ti2T2i^"^ 



To = 



2(1 - T12T21) V12 

2(1 - T2lTi2y^T21 (1 + r2iri2)(l - r2lTi2)"^ 



We the aid of this parametrization, one can show that Q 



(5.27) 



(5.28) 



(5.29) 



where it is understood that a complex conjugate, such as )*, should be expressed in terms of ^21"- It 
follows from this that the result for the invariant integration measure on the GUE coset space in terms 
of eigenvalues and angles is given by 



(Ao - Ai)- 



r c?Ao A 1 d(f>Q d(j)idridrf dpdp* 



(5.30) 



Combining this with Eqs. (4.74), (5.21) and (5.25), we obtain the result 

1 



N-^G{Et)G{E2) = 



4A2 



VQ trgAufcu-HrgAwfcw-ie2*"*'g^ 



/ PQ(Ao-Ai)^*77P>e2-(^''"^^^ 



dAo / dAie2'"(^«-^i) 
1 J -1 
2 



1 

A2 
2i 

A2 \ t,lj+ 



/d 



sm 



■KLJ ' 



(5.31) 



where the abbreviation conn, stands for 'connected part', and recalling that 2a — ttlo^ /d. We should point 
out that the integral over the contribution to Eq. ( ^.31 ) of zcroth order in the Grassmann variables, namely 
(Ao + Ai)2, is singular and actually non- vanishing when treated carefully. However, this contribution 
simply represents the disconnected part of the correlator, which we want to subtract in any case. 



To see this, we appeal to a theorem quoted in Ref. 23 which states that the contribution to the super- 
integral J'DQf{Q) from the component of f{Q) that is of zeroth order in the Grassmann variables 
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is simply given by f{L), i.e. the function f{Q) evaluated at the origin of the Q-coset manifold. The 



corresponding contribution from the integral in Eq. (4.74) is thus seen to be 



-(^trgfcLn)-(^trgfcL22 



(5.32) 



i.e., it is just the disconnected part. Let us write C (a;) = C{Ei, E2)\e^q- Then we obtain, fromEq. ( 2.12D , 
the expression 



M)2 



(5.33) 



for the connected part of the density-density correlator at the centre of the spectrum. The (5-function 



arises from the singular contribution to the real part of Eq. (5.31) when uj 0, which is induced by the 
positive imaginary part of lu^ = lu + ie. The corresponding term is given by 

1 1 



2d I I 



2d 1 



2d 

+0+ A' 



(5.34) 



A description of the modifications required of the present formalism in order to deal with the GOE 
is given in Ref. ^ The presence of time-reversal symmetry (manifested in the real-symmetric nature of 
the Hamiltonian H) necessitates a doubling of the dimension of the supervector ip before performing the 
Hubbard-Stratonovich transformation, viz. <i> = {(p, s(p*)^ , so that a Majorana-type condition $* = C<i> is 
satisfied, where C is an appropriate 'charge- conjugation' matrix. This entails a doubling of the dimension 
of the graded matrices Q, supplemented by a reality condition Q* = CQG^^. 
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APPENDIX A 

In this appendix, we collect the various constant matrices that appear throughout the paper and we 
present derivations of the matrix Gaussian integrals over commuting and anti-commuting variables that 
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appear in Eqs. mM) and (4.7), respectively. Matrices arising from the graded structure of the theory and 



the symmetry-breaking due to the presence of advanced and retarded components are given by 

K - I -1 i ' ~ I -1 i ' '^PP ~ [ ~k ■ ^^■'^> 



pp' 



Matrices associated with the source terms are given by 



V(l)=( J o) - V(2)=(!! ?) . (A.2) 
V / pp' \ / pp' 

In all cases, the explicit indices indicate the space to which the displayed block structure pertains. 



To derive Eq. (4.6), we first observe that, since the matrix A that appears in the exponent is assumed 
to be positive Hermitian, it can be diagonalized according to A = WAbU, where U is a unitary matrix 
and Ad is the diagonal matrix of the positive eigenvalues, Aj^ — diag(Ai, A2, . . . , Ajv). This decomposition 
facilitates the linear transformation of integration variables 

z' ^Uz , z^ ^ z'^ = z^U'^ , (A.3) 

whose Jacobian is equal to unity. The resulting decoupling of the integration variables z' gives rise to a 
factorized form of the Gaussian integral in Eq. ( |4.6| ): 

N 

HA) ^\[^ (dz'J dzi* e-'^^'"'^'" . (A.4) 
fe=i 

Next, we pass to the real and imaginary parts of the holomrphic variables z'j. — Xk + iyk, z'j.* = Xk — iyk, 
in which case 



h{A) = W- dxke-^"^^ / dyke 

I ^ ^ J — QO J —00 



k=l 
N 



fe=i 

= (detA)-i , (A.5) 
having appealed to the fact that the product of eigenvalues coincides with the determinant. 



The analogous integral over anti-commuting variables Eq. (4.7) is most easily derived by expanding 
the integrand in a power series. Since rjf — 77^^ = and we have 27V Grassmann degrees of freedom {rji, rj^ 
for £ = 1,2,..., N), only terms up to order 2N can survive this power-series expansion of the exponential. 
On the other hand, in any of the ensuing polynomial terms that are of order less than 2N , at least one 
Grassmann degree of freedom (say rjk) will be missing. Since, by Eq. ( [4.9| ), J drjk — 0, it follows that such 
terms vanish on integration. Thus, only the polynomial of order 2N survives; and so we can write for 



the Gaussian integral in Eq. (4.7) 



IM)- ifi'^^-^^iv^Ar^- . (A.6) 

=1 
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Now, 

N 

Tj^Arj = J2 "nlA.jT^, . (A.7) 

Again, since r^j = ry^^ = 0, we have 

N 

<T,TePjv i=i 

SO that each rji and appears exactly once for all i = 1, 2, . . . , TV, and where cr, r range over the group 
of permutations of N objects Pat. 

After some elementary manipulation of the product of Grassmann variables above to achieve the 
relation 

N N N 

e=i 1=1 1=1 

where (—1)'^ equals +1 or —1 if the permutation a is even or odd, respectively, we can recast the integral 



(A.6) into the form 

N N 



o-,rePjv 1=1 
N 



^ ^ (-ir(-irn^-w-W' (A.io) 



m 

a-,TePN 1=1 



having made use of Eq. (4.8). Next, let us write £ = a ^(fc) for appropriate k — 1, 2, . . . , A^. Then 



N N N 



Hence, 



^(t(£),t(j') — n ^o-(o-i(fc)),T((j-i(fe)) — JJ^ ^A;,ro<T-i(fc) ■ (A. 11) 

e=l k=l k=l 



1 - ^ 



A! 

o-ePivreFjv k=l 



(k) 

k=l 

m E (-irn^^-M.) 

cr,Z''£Pjv /e— 1 

E (-irn^w 

i/ePiv fe=i 

E (~l)''^l,!^(l)"^2,i/(2) • ■ • ^AT 
i/GPjv 

detA (A. 12) 



by the definition of the determinant. 
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APPENDIX B 



If we write for the block-diagonal matrix R introduced above Eq. (4.48), i? = diag(i?i, i?2), then 
R^^ = sR^s implies that R^^ — r\ and i?2^^ — kR^k. Since Ri is a unitary 2x2 supermatrix, we 
can parametrize it as 

i?i=exp{^(f )}exp( -f^y (B.l) 
Then, since 5P = R~^PdR, we have 

1 - iera d \f a, o \f i- icra -a ^ ^ 

(B.2) 



ai ~ {ai - ibi)^l^i {ai ~ ibi)Q 

(ai - ibi - (ai - 

Since i?2 is a pseudo-unitary 2x2 supermatrix, we can parametrize it as 



It follows that 



R, = exp^ ^ ( j | expjz « ) [ . (g.S) 



^ \^ (62 + «a2)42 «02 + (02 - «02)C242 ' ^ ^ 



The nilpotent parts on the diagonal can be eliminated by appropriate 'contour deformations'. Then 
in both cases p = 1, 2, we have the form 

with 

VPp = -idP^°dP°^dP^°dPl^ 
= dx dy dfj drj 

= —rdrdOdrjdfj. (B.6) 
Thus, trg(5Fp)^ — r'^ ~\- 27777, and so, from Eq. ( 4.69|) , 



2lT fOO 



Kp = - de rdre-^""'' '"^ / di^ df) e-"^^"^ = 1 , (B.7) 



/o ^0 

recalling the normalization of Grassmann integrals in Eq. ([l.§|). 



APPENDIX C 



A direct representation of the density of states can be obtained by taking for the source matrix 



J ^eM ^ek® L. Then, from Eq. ([4.12D , 



= -TrgfcL£)-i 



E = 
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= -2 [Tr G+ [E) - Tr G~ {E)] 
^ AiirNpiE) , 



(C.l) 



noting that here we should take Co = e. From Eqs. (4.71) and ( 4.73| ), we have 



e=0 



— j'DQ trgkLQ exp|-— trgiQj , 



(C.2) 



which imphes 



p(0) = -1- hm / VQ tYg kLQ exp|-^trgLQ| . (C.3) 
Clearly, only the part of the integrand that is of zeroth order in the Grassmann variables contributes. 



Thus, by the theorem cited below Eq. ( 5.31 ), we obtain 

1 



p(0) = — trgfcL2 = — 
47rA ttA 



(C.4) 



Also, from Eq. (4.22) and (4.23), 



1^ 



E = 



Thus, 



N 

<tTgkLS>2 = —<tTgkLQ>3 = AnNp 
A 



(C.5) 



(C.6) 



This establishes the relation between the expectation of Q in the £3-theory and S in the >C2-theory. It is 
clear that <Q>3 = L. 



Another expression for the average density of states can be obtained by appealing to the so-called 
Ward identity, 

1 „ 1 



G+iE)-G-{E) 



Tr- 



Tr- 



E-H + ie E -H-ie 
-2ieTr- ^ 



{E - HY + £2 



in which case 



p{E) — lim e Tr 



1 



ttN e^o+ E~ H + ie E - H - ie 
lim eTr£)r/D„^ . 



(C.7) 



(C.8) 



Here, we simply take uj = e. A source matrix that will generate this combination is furnished by Eq. ( 4.15 ) 
with Mjn = km ® m ^ B.F, where 



. , ,01 



Uaa' (1 

' - 1 

pp' \ / aa' 



uaa _ ( 
' '^F - I 1 



(C.9) 
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Then, 



92 



dSBdSB 



Z{e) 



(C.IO) 



e=0 



A similar source matrix was used in Rcf. |2| for calculating moments of conductance. In the present case, 
we have trgLJ = but trg = 2trg(e-/c)2 = 2{e% ~ e],). Also, 



trgQJ = trg(e-fc)Qi2 +trg(e-/c)Qi2 

= es trgfcB((5i2 + Q21) + e_F trgfcF(Qi2 + Q21) 



It follows that 



whence 



del' 



e=0 



^[trg/cB(Ql2+Q2l)]' 



p(0) ^ - 



TV 
2^ 



^lun e j VQ [trg fc^ (Q12 + Q21)]' exp{-^ trg Lq} . 



(c.ii) 



(C.12) 



(C.13) 



We see immediately that only the term in the integrand of highest order in Grassmann variables con- 



tributes to this expression. Since, from Eqs. (5.7) and (5 



we have 



Q12 = 2u ^^y/T+Jlflv , Q21 — 2v ^^iy/l + fifiU 



trgfci3(Qi2 + Q2i) = 2[(^o+Mo)(l + Mo/Io)(l-5'7*^)(l + l//') 



(C.14) 



Thus, the contributing part of the integrand is given by 



(C.15) 



[trg/cB(Qi2 + Q2i\ 



l2 c.p. 



-8[(1 + AioMo)AioMo - (1 + MiMi)miMi 
f ^{1 + no'Po){fJ-l+'Po)]v*VP*P ■ 



Now, 



(1 + AioAio)MoMo - (1 + = i(Ao - A^) 



(C.16) 



(C.17) 



while the third term in the square brackets above produces a vanishing contribution after integration 
over the angles (j)o,(j)i. Hence, we obtain 

^ = -i- hm e' r dXo r dX, ^i±^e-'(^-^^) , (C.18) 
SttA £'-*o+ Ji J_i Ao - Ai 

on recalling that d — nX/N and setting e' — 2'Ke/d. We can write the integral in Eq. ( |C.18 ) as 

dz J dXo J rfAi (Ao + Ai)e-^(^°-^^^ (C.19) 

The variable transformation 



Ao = 1 + 2to , Ai = 1 - 2ti 



(C.20) 
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yields the form 



OO pi 



dz / dtie-^'*' / dto[(l-ii)+to]e"^^*" 
Jo Jo 



2 



(C.21) 



This leads immediately to the expected result p(0) — l/(7rA). 



Let us explain why we refer to Eq. ( p. 71) as a Ward identity. With a slight modification of our 



notation, we have from Eqs. (4.71) and (4.73) 



Z{J + ii:}L)^ JvQ exp|^trgg(J + iwi)| . 



Thus, we can also write 



Z{A) = Jt^Q expj^trgQ^} 
for any supermatrix A. Now let A = Aq + 5A where 5A is infinitesimal. Then we have 

7 AT r riN -i 

Z{Ao + 6A)~Z{Ao)^— / expj— trgQ^I trgQSA . 

Our theory has the symmetry property 



z(r-i( J + ic;)L)T) = z{j + iCjL) , 



(C.22) 



(C.23) 



(C.24) 



(C.25) 



which is evident both from Eq. (4.22) and Eq. (4.71). For infinitesimal transformations, T — 1 + dT^ 
= 1 - this reads 



Z{J + iuL +[J + iCoL, 5T]) - Z{J + iQL) = . 



So if we choose = + ii^i, 5 A = [A^, ST], then we obtain the Ward identity 



noting that 

trgQ^yl 

and 6T is arbitrary. Equivalently, we can write 

iN 



trgST[Ao,Q] 



iN 



(C.26) 



(C.27) 



(C.28) 



(C.29) 



We can use Eq. (C.29) to generate families of integral identities. For example, since on the RHS 

Q 



(C.30) 
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we have 



trgkc,ML[L,Q] ^ 2trgfc„(gi2 + Q21) 
tTgkaML[J,Q] = 2tTg[kaML,J]Q 
= 2trg 



[kaJu + J2lka) 
kaJll — J22ka 



-kaJ22 + Jllka 
kaJl2 + J2lka 



for a = B,F, etc. If we suppose that [ka, Jpp'] = 0, then 



trgkaML[J,Q] = -trgfc„(Ji2 + J2i)£(9 + trg fc„( Jn - J22)MQ , 



and the Ward identity (C.29) now reads 

j VQ [trgA:,(Jn- J22)(Qi2 + 02i)-trgA;„(Ji2 + J2i)(Qii-Q22)] eT^ 



iN 

-2iCj I VQ trgfc„(Qi2 + Q2i)e-*'-*5'5(^+'"'^) 



Let us also take J = epkpL + e'^k^Ai, so that 

Jll - J22 = 'isfikfj , J12 + J21 = ^e'^k^ , 

while 

trgQJ = eptrgki3{Qii - Q22) + e'^ trgk^{Qi2 + Q21) ■ 
Then setting a = /3 = 7 and equating powers of e and e' yields the family of relations 

F{m + 1, n + 1) = (m + l)F{m, n + I) - nF{m + 2, n - 1) , 

valid for n = 0, 1, 2, . . . and m = —1, 0, 1,2,.. ., where 

F{m,n)^ / VQe-—''^'^'^ [trgfc„(Qn - Q22)]" [trgfcc(Qi2 + Q21)]" 



(C.31) 



(C.32) 



(C.33) 



(C.34) 



(C.35) 



(C.36) 



(C.37) 



This constitutes the basis for a recurrence relation for the Q-integrals F(m, n). The case m = — 1, n = 1 



yields an identity equivalent to Eq. (C/7). We see that the Ward identity relates the transverse modes 
Q12, Q21 to the longitudinal modes Qii, Q22- 
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